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We review the direct detection of supersymmetric dark matter in the light of recent 
experimental results. In particular, we show that regions in the parameter space of several 
supergravity scenarios with a neutralino-nucleon cross section of the order of 10 -6 pb, i.e., 
where current dark matter detectors are sensitive, can be obtained. These are scenarios with 
large tan /3, with non-universal soft super symmetry-breaking terms, with multi-TeV masses 
for scalar superpartners known as 'focus point' supersymmetry, and finally scenarios with 
intermediate unification scale which appear naturally in some superstring constructions. 

1 Introduction 



One of the strongest motivations for physics beyond the standard model is the exis- 
tence of dark matter. Substantial evidences exist suggesting that most of the mass 
in the Universe is some non-luminous matter, the so called 'dark matter', of as yet 
unknown composition. Currently the most convincing observational evidence for the 
existence of dark matter comes from the analysis of rotation curves of spiral Galaxies, 
i.e., measurements of the velocity of isolated stars or gas clouds which orbit in the 
outer parts of spiral Galaxies. It has been noted that there is not enough luminous 
matter in those Galaxies to account for their observed rotation curves [l[|. 

The detailed analysis of rotation curves offers convincing evidence that 90% or 
more of the mass of the Galaxies is dark. Other evidences arise from large scale 
measurements, as e.g. the motions of cluster members Galaxies. A key question is 

*Extension of a talk given at Cairo International Conference on High Energy Physics (CICHEP 
2001), January 9-14. 
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then, what could this dark matter be? We don't know the answer yet, but we do know 
that not all of it can be ordinary (baryonic) matter, since measured abundances of 
helium, deuterium and lithium in the scenario of big-bang nucleosynthesis impose a 
strong upper bound on the baryon density in the Universe. This density is too small 
to account for the whole dark matter in the Universe. The conclusion is that baryonic 
objects, such as e.g. MACHOs, can be components of the dark matter, but more 
candidates are needed ||. 

Fortunately, particle physics offers various candidates for (non-baryonic) dark mat- 
ter, all of which would indicate new physics beyond the well tested standard model 
of particle physics [[|. For example, long-lived or stable weakly-interacting massive 
particles (WIMPs) can remain from the earliest moments of the Universe in sufficient 
number to account for a significant fraction of relic density. These particles would form 
not only a background density in the Universe, but also would cluster gravitationally 
with ordinary stars in the galactic halos. 

This raises the hope of detecting relic WIMPs directly, by observing their elastic 
scattering on target nuclei through nuclear recoils. Since WIMPs interact with ordinary 
matter with very roughly weak strength, and assuming that their masses are of the 
order of weak scale (i.e., between 10 GeV and a few TeV), it is natural to expect a 
WIMP-nucleus cross section of the same order as that of a weak process, which is around 
1 pb. This would imply a WIMP-nucleon cross section around 10~ 8 pb, too low to be 
detected by current dark matter experiments, DAMA, CDMS and UKCDM, which are 
sensitive to a cross section around 1CT 6 pb. Surprisingly, the DAMA collaboration 
reported recently 0] data favouring the existence of a WIMP signal in their search for 
annual modulation. When uncertainties as e.g. the WIMP velocity or possible bulk 
halo rotation, are included, it was claimed that the preferred range of parameters is (at 
4a C.L.) KT 6 pb < a < 10" 5 pb for a WIMP mass 30 GeV < m < 200 GeV. Unlike 
this spectacular result, the CDMS collaboration claims to have excluded || regions of 
the DAMA parameter space. 

In any case, due to these and other projected experiments, it seems very plausible 
that the dark matter will be found in the near future. In this situation, and assuming 
that the dark matter is a WIMP, it is natural to wonder how big the cross section 
for its direct detection can be. The answer to this question depends on the particular 
WIMP considered. The leading candidate in this class is the lightest neutralino ||, a 
particle predicted by the supersymmetric (SUSY) extension of the standard model. In 
this paper we critically reappraise the known SUSY scenarios based on neutralinos as 
dark-matter candidates, and in particular the scenarios constructed recently in order 
to enhance the neutralino-nucleon cross section. This is the case of scenarios with large 
tan/5 with non-universal soft SUSY-breaking terms [0, j|, [T0| , with multi-TeV 

masses for scalar superpartners known as 'focus point' supersymmetry [[□]], and finally 



scenarios with intermediate unification scale |12|-p4 
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2 Supersymmetric predictions for the neutralino-nucleon cross 
section 



In SUSY models, -R-parity is often imposed to avoid weak scale proton decay or lepton 
number violation. Imposing this symmetry yields remarkable phenomenological impli- 
cations. SUSY particles are produced or destroyed only in pairs and, as a consequence, 
the lightest supersymmetric particle (LSP) is absolutely stable. The former implies 
that a major signature for .R-parity conserving models is represented by events with 
missing energy (for instance, e + e~ — > jet + missing energy). The latter implies that 
the LSP might constitute a possible candidate for dark matter. Concerning this point, 
it is remarkable that in most of the parameter space of SUSY models the LSP is an 
electrically neutral (also with no strong interactions) particle, called neutralino. This 
is welcome since otherwise the LSP would bind to nuclei and would be excluded as a 
candidate for dark matter from unsuccessful searches for exotic heavy isotopes ]15 . 



In the simplest SUSY extension of the standard model, the so-called minimal super- 
symmetric standard model (MSSM) there are four neutralinos, xl = 1)2,3,4), since 
they are the physical superpositions of the fermionic partners of the neutral electroweak 
gauge bosons, called bino (B°) and wino (W®), and of the fermionic partners of the 
neutral Higgs bosons, called Higgsinos (H®, H^). Therefore the lightest neutralino, Xii 
will be the dark matter candidate. The neutralino mass matrix with the conventions 
for gaugino and Higgsino masses in the Lagrangian, C — \ YJ a M a X a X a + [iH®H® + 
h.c, is given by 



(Mi -M z cos/3sm9 w M z sm.(3sm.0 w \ 

M2 Mz cos (3 cos 6w — M z s\n (3 cos 9w 

— Mz cos f3 s\n 9 w Mz cos (3 cosOw — fi 

y Mz sin (3 sin Ow — Mzsm. (3 cos 9w — M j 



(1) 



in the above basis (B° = = —iX 3 , H®, H%). Here M x and M 2 are the soft bino 

and wino masses respectively, /i is the Higgsino mass parameter and tan f3 = (H®) / (H^) 
is the ratio of Higgs vacuum expectation values. We parameterize the gaugino and 
Higgsino content of the lightest neutralino according to 

X° = N n B° + iV 12 H>3° + N 13 H° d + N U H° U . (2) 

It is commonly defined that Xi is mostly gaugino-like if P = |iVn| 2 + |iV 12 | 2 > 0.9, 
Higgsino-like if P < 0.1, and mixed otherwise. 

The relevant effective Lagrangian describing the elastic Xi-nucleon scattering in the 
MSSM is given by 

£e// = { a gXXQQ + (3 q xi^ 5 xqi^ 5( i) > ( 3 ) 

q 

where a q , f3 q are given e.g. in ref.Jl6| with the sign conventions for Yukawa couplings 



in the Lagrangian, C = —h u H^ULUR — hdH^didu — h e H®eLen + h.c, and the sum runs 
over the six quarks. The contribution of the scalar (spin- independent) interaction, the 
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Figure 1: Feynman diagrams contributing to neutralino-nucleon cross section. 



one proportional to a q , to the x'J-nucleon cross section is generically larger than the 
spin-dependent interaction, the one proportional to f3 q . One can concentrate then on 
the scalar x^-nudeon cross section. This is given by 



a = , (4) 

7T 



where m r is the reduced Xi mass and 



; £ /.A (5) 



m g 'm q 



with m the mass of the nucleon and m q the mass of the quarks. The scalar coefficients 
a q include contributions from squark (q) exchange and neutral Higgs (h, H) exchange, 
as illustrated in Fig. [I]. The numerical values of the hadronic matrix elements f q for 
the proton are as follows | 16| : 



f u = 0.020 ± 0.004 , f d = 0.026 ± 0.005 , f s = 0.118 ± 0.062 , (6) 

and f c ^t = §f(l — J2 q =u,d,s fq) ■ F° r the neutron the value of f s is the same and 
therefore the scalar cross sections for both, protons and neutrons, are basically equal. 

The cross section for the elastic scattering of relic neutralinos on protons and neu- 
trons has been examined exhaustively in the literature |J. This is for example the 
case in the framework of minimal supergravity (mSUGRA). Let us recall that in this 
framework one makes several assumptions. In particular, the scalar mass parameters, 
the gaugino mass parameters, and the trilinear couplings, which are generated once 
SUSY is broken through gravitational interactions, are universal at the grand unifica- 
tion scale, M G ut ~ 2 x 10 16 GeV. They are denoted by m , Ml/2, and A respectively. 
Likewise, radiative electroweak symmetry breaking is imposed, i.e., \fi\ is determined 
by the minimization of the Higgs effective potential. This implies 



rri 2 Hd - m 2 Hu tan 2 g _ 1 
tan 2 (3-1 2 



fj, = T §^- — - - -M| . (7) 



With these assumptions, the mSUGRA framework allows four free parameters: mo 
Mi/2, A , and tan/?. In addition, the sign of \x remains also undetermined. 

a Larger values, as for example f s = 0.455, have also been used in the literature, see e.g. ref. B. 
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Figure 2: Running of the soft Higgs masses- squared with energy. 

It was observed (for a recent re-evaluation see ref. |16|]) that for low and moderate 
values of tan/5 the predicted scalar neutralino-proton cross sections are well below 
the accessible experimental regions. To understand this result qualitatively, we show 
schematically in Fig. |2] the well known evolution of m 2 H and rn 2 Hu with the scalef]. 
Since m 2 Hu becomes large and negative (notice also that m 2 H , neglecting bottom and 
tau Yukawa couplings, becomes positive), due to relation (□) becomes also large, 
in particular much larger than M\ and M 2 . Thus, as can be easily understood from 
eqs.(|I]) and @, the lightest neutralino will be mainly gaugino, and in particular bino 
since at low energy Mi « \M 2 . 

We show this fact in Fig. ||a, where for tan (3 = 3 the gaugino- Higgsino components- 
squared iV^ of the lightest neutralino as a function of its mass m^o are exhibited. Here 
we are using as an example mo = 150 GeV, and M\ii is essentially fixed for a given 
m^o. Note that Nx% is extremely large and therefore P > 0.9. Then, the scattering 
channels through Higgs exchange shown in Fig. [I] are suppressed (recall that the Higgs- 
neutralino-neutralino couplings are proportional to -/V 13 and N u ) and therefore the cross 
section is small. As a matter of fact, the scattering channels through squark exchange, 
shown also in Fig. [I], are also suppressed by the mass of the squarks. Indeed in this 
limit the cross section (|4]) can be approximated as 

a 2 

cr^o-p ~ m r — — ^ f s \N U 



Thus it can be roughly estimated to be cr$o_ p ~ 10 -8 pb, for a ~ 10 -2 , f s ~ 10 _1 , 
m r ~ 1 GeV and m s ~ 300 GeV. This is precisely the value of the cross section 
expected for a weak process as discussed in the Introduction. 

In Fig. |3|b we show the cross section function of the neutralino mass m^o 

for tan/5 = 3. We also choose to plot only the case \i > since for negative values of n 
the cross sections are much smaller[|. Fig. [3|b has been obtained using formula (|j) with 

b Let us remark that, given the convention used throughout this review for gaugino masses in the 
Lagrangian, C = ^ T\ Mg^a^a + h- c -> we are using the renormalization group equations (RGEs) 
obtained e.g. in ref. |17| , but with an opposite sign in the gaugino contributions to the RGE's of the 
Aq parameters. 

c It is also worth noticing that constraints coming from the b — > sj process highly reduce the /i < 
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Figure 3: (a) Gaugino-Higgsino components-squared of the lightest neutralino as a 
function of its mass, (b) Scatter plot of the neutralino-proton cross section as a function 
of the neutralino mass for the scenario discussed in the text. DAMA and CDMS current 
experimental limits and projected GENIUS limits are shown. 

the central values for the hadronic matrix elements given by eq.(^). We will use these 
values throughout this review^. Likewise we impose the present bounds coming from 
accelerators. These are LEP and Tevatron bounds on SUSY masses and CLEO b — > S7 
branching ratio measurements. For the Higgs mass we use the present experimental 
lower limit for tan/3 = 3 Hl~9"l , i.e. rrih > 95 GeV. Although, for the values of tan/? that 



we will consider in the next section, tan/5 = 10, this limit is considerably weaker ]19 
rrih > 80 GeV, we will still use > 95 GeV. Concerning the parameter space of the 
figure, we have used 30 GeV < mo < 550 GeV, where the lower bound is in order to 
avoid the stau being the LSP for all values of m^o. Note that the curve associated to 
tuq = 550 GeV corresponds to the minimum values of cr^o_ p in the figure. On the other 
hand, since the cross sections are not very sensitive to the specific values of A in a 
wide range (we have checked that this is so for | A /Mi/ 2 | < 1), we fix A = Mx/ 2 in 
the figure. As mentioned above the gaugino mass M1/2 is essentially fixed for a given 
tan/3 and m^o. For this figure we have to take 140 GeV < Mi/ 2 < 350 GeV, where 
the lower bound, corresponding to m^o > 55 GeV, is due to the experimental bound 
on the lightest chargino mass m~± > 90 GeV. Negative values of Mi/ 2 , not shown in 
the figure, correspond to smaller cross sections. 

As we can see in the figure cr^o-p ~ 10~ 7 pb, an d therefore for these values of the 



parameters we would have to wait in principle for the projected GENIUS detector |20 



parameter space. In addition, recent measurements of the anomalous magnetic moment of the muon 
[ PI seems to exclude all cases with fx < 0. 

d It is worth mentioning here that values for the hadronic matrix elements as those discussed in 
footnote a will increase in one order of magnitude the results of this and the other figures in this 
review, since their contribution appears raised to the square in eq.fl). 
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Figure 4: The same as in Fig. |3]b but for tan/3 = 10 and 20. 

to be able to test the neutralino as a dark-matter candidate. 

3 Large neutralino-nucleon cross section in supergravity sce- 
narios 

Recently there has been some theoretical activity [f7|-|14| trying to obtain regions in the 
parameter space of supergravity (SUGRA) scenarios compatible with the sensitivity 
of current dark matter detectors, DAMA and CDMS. The key point in most of these 
scenarios to carry it out consists of reducing the value of Following the discussion 
below eq.([?D in the previous section, the smaller is, the larger the Higgsino com- 
ponents of the lightest neutralino become. Eventually, will be of the order of Mi, 
M.2 and Xi will be a mixed Higgsino-gaugino state. Indeed scattering channels through 
Higgs exchange in Fig. [TJ are now important and their contributions to the cross section 
@ can be schematically approximated as 

^ ? _ p ^m^/ s 2 |iV 1 3iV 11 | 2 , (9) 
m h 

where X s is the strange quark Yukawa coupling and mu represent the Higgs masses. 
With the same rough estimate as in eq.®, one obtains cr^o_ p ~ 10~ 6 pb for m/, ~ 100 
GeV. 

In this section we review the different SUGRA scenarios that can be found in the 
literature in order to enhance the neutralino-proton cross section cr^o_ p to be of the 
order of 10 -6 pb. 
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3.1 Scenario with large tan/3 

In the previous section we showed in Fig. ^|b the neutralino-proton cross section for 
tan/3 = 3, in the framework of mSUGRA. Here we show the same in Fig. |] but for 
larger values of tan/3, in particular for tan/3 = 10 (green points) and 20 (red points). 
We see that the cross section increases when the value of tan /3 increases. For moderate 
values of tan/3, the reason being that the top(bottom) Yukawa coupling which appears 
in the RGE for iTi"H u {m 2 Hd ) decreases(increases) since it is proportional to ^~^(~g)- 
This implies that the negative(positive) contribution m 2 H {m 2 H ) to fi 2 in eq.(|7|) is less 
important, and therefore decreases. 

However, for tan/3 > 10, the value of fi 2 is very stable with respect to variations of 
tan f3. This is due to the fact that /z 2 ~ —m 2 H2 — \M\ (see eq.(|?])). Since sin f3 ~ 1, the 
top Yukawa coupling is stable and therefore the same conclusion is obtained for m 2 H2 
and /i 2 . Thus, the reason for the cross section to increase when tan/3 increases cannot 
be now the increment of the Higgsino components of the LSP. Nevertheless there is a 
second effect in the cross section which is now the dominant one: the contribution of 
the down- type quark Yukawa couplings (see eq.([])) which are proportional to 

It was in fact pointed out in refs.|7|, [8| that the large tan/3 regime allows regions 
where (J^o_ p « 10~ 6 pb is reached. We can see in Fig. f| that this is so for tan /3 > 20. 

Very large values of tan /3, like tan /3 ~ 50 which correspond approximately to the 
unification of the tau and top Yukawa couplings at M GUT , have also been considered ||. 
Although it was found, as expected, that the cross section is enhanced, the well known 
experimental limits coming from b — > 57 for large tan /3, lead to severe constraints on 
the parameter space. In particular, these constraints imply cr^-p ~ 10~ 6 pb. 



3.2 Focus point supersymmetry scenario 

Another possibility to obtain a large neutralino-nucleon cross section arises in the so- 
called focus point supersymmetry scenario. This has been proposed |2T| in order to 
avoid dangerous SUSY contributions to flavour and CP violating effects. The idea is to 
assume the existence of squark and sleptons with masses which can be taken well above[] 
1 TeV. It has also been argued that this situation produces no loss of naturalness. 
The implications of focus point supersymmetry for neutralino dark matter have 



been considered in ref.|Tl[]. In particular, it was pointed out that for m > 1 TeV, 
unlike the usual cases with m < 1 TeV, the lightest neutralino is a gaugino-Higgsino 
mixture over much of parameter space. Let us recall from the previous subsection that 
for moderate and large values of tan /3, /x 2 in eq. (^) can be approximated as 

/i 2 « -m 2 Hu - l -M 2 z . (10) 

Thus for m > 1 TeV, m 2 H becomes less negative, and therefore decreases. (For 
tan/3 < 5, /i becomes sensitive to ^# d , and as m increases also increases, and 

e Of course this scenario rules out SUSY as an explanation of the possible deviation in the muon 
anomalous magnetic moment from the standard model prediction |18| . 
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Figure 5: The same as in Fig. but for tan/? = 10 and using non-universal soft 
scalar masses as explained in the text. The universal case (green region) is shown for 
comparison. 



so there is no mixed gaugino-Higgsino region for small tan/3.) However, as discussed 
in ref. ITT], one still needs large values for tan/? (of the order of 50) in order to have 
pb. 



3.3 Scenario with non-universal soft terms 

The soft SUSY-breaking terms can have in general a non-universal structure in the 
MSSM. Such structure can be derived from SUGRA p^[ . For the case of the observable 
scalar masses, this is due to the non-universal couplings in the Kahler potential between 
the hidden sector fields breaking SUSY and the observable sector fields. For the case 
of the gaugino masses, this is due to the non-universality of the gauge kinetic functions 
associated to the different gauge groups. Explicit string constructions, whose low- 
energy limit is SUGRA, exhibit these properties . 

It was shown recently, in the context of SUGRA, that non-universality allows to 
| [10]] the neutralino-proton cross section for moderate values of tan (3. This 



increase 



can be carried out with non-universal scalar masses and/or gaugino masses, as we will 
discuss below. 



(i) Non-universal scalar masses 



Non-universality in the Higgs sector, concerning dark matter, was studied in refs. [ 
24| . Subsequently, non-universality in the sfermion sector was added in the analysis 
PS] , [5]. In order to avoid potential problems with flavour changing neutral currents, 
it was assumed that the first two generations of squarks and sleptons have a common 
scalar mass m at Mqut-, and non-universalities were allowed in the third generation 
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and Higgs sector. Thus the soft masses can be parameterized as follows: 



m H d = "io( 1 + 5 i)> m H u =mjj(l + 5 2 ), 

m Q L = m 2 (l + 6 3 ), m 2 R = mjj(l + <? 4 ), 

m L = ^(l + 5 5 ), m d 2 R = m 2 {l + 8 e ), 

m l L = rn 2 (l + 5 7 ), 



(11) 



where Ql = (t L ,b L ), L L = (vl, t l), u r = t R , e R = t r , and it is assumed that —1 < 
Si < 1. 

As explained at the beginning of this section 3, fi 2 is one of the important param- 
eters when computing the neutralino-nucleon cross section. Its value is determined by 
condition (0) and can be significantly reduced for some choices of 5's. We can have 
a qualitative understanding of the effects of the 5's on \i from the following. First, 
when m 2 Hu (m 2 Hd ) at M GUT increases(decreases) its negative (positive) contribution at 
low energy in eq.(|7|) is less important. Second, when m 2 Q L and m 2 R at M GUT decrease, 
due to their contribution proportional to the top Yukawa coupling in the RGE of m# 
the negative contribution of the latter to fi 2 is again less important. Thus one can 
deduce that fi 2 will be reduced (and hence c^,o_ p increased) by choosing <5 3 , 5 4 , 5\ < 
and <5 2 > 0. 

Following this analysis [0, [£| we show in Fig. [5|, for tan/3 = 10, a scatter plot of 

(X^o p as a function of m^o for a scanning of the parameters as follows: < 8% < 1, 

— 1 < ^1,3,4 ^ 0- The other 5's are not so important in the computation and we take 
them vanishing. We are also taking A = M ly / 2 , 30 GeV < tuq < 550 GeV as in Fig. |]b, 
and we have to take now 140 GeV < Mi/ 2 < 450 GeV. For comparison we superimpose 
also the region (green area) obtained in Fig. |] for tan/? = 10 with universality, 5i = 0. 
We see that non- universal scalar masses can help in increasing the values of <J^o_ p . In 
fact non- universality in the Higgs sector gives the most important effect, and including 
the one in the sfermion sector the cross section only increases slightly. Of course, as 
discussed in detail in Subsection 3.1, for larger values of tan/3 one can get larger cross 
sections. 

(ii) Non-universal gaugino masses 

The effects of the non-universality of gaugino masses on dark matter in SUGRA 
scenarios have been studied in ref. PBl [LT], |27fl . In particular, in ref.[|I(J, the authors 
analyze in detail the neutralino-proton cross section in the case of SU(5) unified models 
with interesting results, finding that the cross sections can increase. 

Here we are interested in the case of the standard model gauge group in order to 
compare our results with both, the universal case discussed in Sections 2 and 3.1 and 
the case with non-universal scalar masses discussed above. Let us then parameterize 
the soft gaugino masses at Mqut as follows: 



where M 12) 3 are the bino, wino and gluino masses, respectively, and we assume that 



Mi 



1 



M 1/2 (l + 5[), M 2 = M 1/2 (l + 5' 2 ), M 3 = M 1/2 (l + 5' 3 ) 



(12) 



-1<*{<1. 
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Figure 6: The same as in Fig. |3|b but only for tan /3 — 10 and using non-universal soft 
gaugino masses as explained in the text. The universal case (green region) is shown 
for comparison. 



Let us discuss now which values of the parameters are interesting in order to increase 
the cross section, with respect to the universal case 5^ = 0. In this sense, it is worth 
noticing that M 3 appears in the RGEs of squark masses, so e.g. their contribution 
proportional to the top Yukawa coupling in the RGE of mfj will do this less negative 
if M3 is small. Therefore /z 2 becomes smaller. 

Taking into account this effect, we show in Fig. || for tan/3 = 10, a scatter plot of 

p as a function of m^o for the following scanning of the parameters: — 1 < 5[ 2 < 1 

and —2/3 < 5' 3 < (where the lower bound is due to the experimental bound on gluino 
masses), imposing M 3 < M 1 ,M 2 . We fix My 2 = 225 GeV and A = M x . As in the 
previous figures we are taking 30 GeV < mo < 550 GeV. We also superimpose for 
comparison the region (green area) obtained in Fig. [| for tan f3 = 10 with universality, 
5[ = 0. Clearly non-universal gaugino masses produce a large increment in the cross 
section. 

Let us finally remark that the upper and right hand side regions with large cross 
sections in the figure correspond to points of the parameter space fulfilling M 2 < M\ 
(and M3 close to its lower bound). As discussed in Section 2, for low or moderate 
tan/3 the lightest neutralino is mainly gaugino P > 0.9, and in particular bino since 
Nn ^> Ni 2 . However, relaxing the universality condition for gaugino masses, M 2 < Mi 
will produce Nn < iV 12 , and therefore the lightest neutralino will have an important 
wino component]]. Since the latter couples through SU(2) interactions, unlike the bino 
component which couples through U(l) ones, the cross section increases. 

(Hi) Non-universal gaugino and scalar masses 

f It is worth noticing that this situation also arises in anomaly- mediated SUSY-breaking scenarios, 
where the cross section can also be increased pa, |29( . 



11 



Given the above situation concerning the enhancement of neutralino-proton cross 
sections through non-universality for moderate tan/5, it is worth analyzing in prin- 
ciple the combination of both possibilities, non-universal gaugino and scalar masses. 
However the consequence for the cross section of combining the parameters studied in 
(i) and (ii) is not very different from the one obtained in Fig. ^. We have explicitly 
checked that the qualitative pattern is basically the same. 



3.4 Scenario with intermediate unification scale 

In the above subsections the analyses were performed assuming the unification scale 
Mqut ~ 10 16 GeV, as is usually done in the SUSY literature. However, it was recently 
realized that the string scale may be anywhere between the weak and the Plank scale 



3"0|-||36||. To use the value of the initial scale, say Mj, as a free parameter for the 
running of the soft terms is particularly interesting since there are several arguments 
in favour of SUSY scenarios with scales Mj « io 10-14 GeV. 

First, these scales were suggested |33], to explain many experimental observations 
as neutrino masses or the scale for axion physics. Second, with the string scale of the 
order of 10 10 ~ 12 GeV one is able to attack the hierarchy problem of unified theories 



without invoking any hierarchically suppressed non-perturbative effect [[33], [3JJ] . Third, 
for intermediate scale scenarios charge and color breaking constraints, which are very 
strong with the usual scale Mqut |37|, become less important [[38]. There are other 



arguments in favour of scenarios with initial scales Mj smaller than Mqut- For example 
these scales might also explain the observed ultra-high energy (« 10 20 eV) cosmic rays 
as products of long-lived massive string mode decays J33], [3j] (see ref.[]39| for more 



details about this possibility). Besides, several models of chaotic inflation favour also 
these scales 



Inspired by these scenarios, it was pointed out [|T^]-[]T4| that the neutralino-proton 
cross section cr^o_ p is very sensitive to the variation of the initial scale for the running 
of the soft terms. In particular, intermediate unification scales were considered. For 
instance, by taking Mj = 10 10-12 GeV rather than Mqut, regions in the parameter 
space of mSUGRA have been found where cr^o_ p is in the expected range of sensitivity 
of present detectors, and this even for moderate values of tan/? (tan/3 > 3). This 



analysis was performed [12] in the universal scenario for the soft terms. In contrast, in 
the usual case with initial scale at Mqut, this large cross section is achieved only for 
tan/3 > 20, as discussed in Subsection 3.1. 

Before trying to understand this result, let us discuss what we mean by an interme- 
diate unification scale. Concerning this point two possible scenarios are schematically 
shown in Fig. |7| for the example Mj = 10 11 GeV. In scenario (a) the gauge couplings 
are non universal, oti ^ a, and their values depend on the initial scale Mj chosen. An 
interesting proposal in order to obtain this scenario in the context of type I string mod- 
els is the following. If the standard model comes from the same collection of D-branes, 
stringy corrections might change the boundary conditions at the string scale Mj to 



mimic the effect of field theoretical logarithmic running [41, 42]. Another possibility 
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Figure 7: Running of the gauge couplings with energy, shown with solid lines, assuming 
(a) non universality and (b) universality of couplings at the initial scale Mj. For 
comparison the usual running of the MSSM couplings is also shown with dashed lines. 



giving rise to a similar result may arise when the gauge groups came from different 
types of D-branes. Since different D-branes have associated different couplings, this 
implies the non universality of the gauge couplings. We will discuss this possibility in 
some detail below. 

On the other hand, scenario (b) with gauge coupling unification at Mj, an = a, 
can be obtained with the addition of extra fields in the massless spectrum. For the 
example of the figure these are doublets and singlets under the standard model gauge 
group. 

As we will see below, the values of the gauge coupling constants at the intermediate 
scale will be important in the computation of the cross section, and scenario (a) will 
be more interesting than (b). 

Let us now come back to the issue of the variation of the cross section with the 
initial scale. The fact that smaller initial scales imply a larger neutralino-proton cross 
section can be understood from the variation in the value of [i with Mj. One observes 
that, for tan/? fixed, the smaller the initial scale for the running is, the smaller the 
numerator in the first piece of eq.(^) becomes. This can be understood qualitatively 
from Fig. |2|. Clearly, the smaller the initial scale is, the shorter the running becomes. As 
a consequence, also the less important the positive(negative) contribution rn 2 H (rn 2 Hv ) 
to \i in eq.([7D becomes. Thus decreases. 

As discussed at the beginning of this Section 3, now the Higgsino components of the 
lightest neutralino, A13 and N\± in eq.(§), increase and therefore the spin independent 
cross section also increases [ I2| . This is shown in Figs. |8] and |9|. In fact these figures 



correspond to the scenario in Fig. [7|a with non-universal gauge couplings. For the 
scenario in Fig. 0b, only with tan f3 > 20 one obtains regions consistent with DAMA 
limits. One of the reasons being that now ^(Mj) and ai(Mj) are bigger than in 
scenario (a) and therefore the low-energy bino and wino masses appearing in eq.(|l]) are 
smaller. As a consequence the increment of the cross section is less important. 

It is also worth noticing that, for any fixed value of Mi, the larger tan/3 is, the 
larger the Higgsino contributions become. The discussion concerning this point in 
Subsection 3.1 for Mqut is also valid for any initial scale Mj. 
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Figure 8: Gaugino-Higgsino components-squared of the lightest neutralino as a function 
of its mass for the unification scale, Mj = 10 16 GeV, and for the intermediate scale, 
Mi = 10 11 GeV. 

In Fig. [8], for tan/? = 10 and mo = 150 GeV, we exhibit the gaugino-Higgsino 
components-squared Nf { of the lightest neutralino as a function of its mass m^o for 
two different values of the initial scale, Mj = 10 16 GeV « M GUT and Mj = 10 11 
GeV. Clearly, the smaller the scale is, the larger the Higgsino components become. 
For Mj = 10 11 GeV, e.g. the Higgsino contribution N13 becomes important and even 
dominant for m^o < 140 GeV. 

The consequence of the above results on the cross section is shown in Fig. [], where 
the cross section as a function of the lightest neutralino mass m^o is plotted. In 
particular we are comparing the result for the scale Mj = M GUT studied in Figs. || 
and |] with the result for the intermediate scale Mi = 10 11 GeV. For instance, when 
m^o = 100 GeV, <y^o_ p for M/ = 10 11 GeV is two orders of magnitude larger than for 
Mqut- In particular, for tan/3 = 3, one finds cr^-p ~ 10~ 7 pb if the initial scale is 
Mj = 10 16 GeV. However cr^o_ p < 10~ 6 GeV is possible if Mi decreases. 

As mentioned before, the larger tan/3 is, the larger the Higgsino contributions 
become, and therefore the cross section increases. For tan/3 = 10 we see in Fig. |9| that 
the range 70 GeV < m^o < 100 GeV is now consistent with DAMA limits. 



Let us remark that these figures have been obtained |C| taking 30 < mo < 550 GeV 
and A = M 1 / 2 as in previous sections. In any case, as mentioned also in Section 2, 
the cross section is not very sensitive to the specific values of A Q . In particular it was 
checked that this is so for | Aq/Mi/2 \ < 1. For example, relation Aq = —M\/2 is 



particularly interesting since it arises naturally in several string models [^, PH . 

Let us finally recall that the above computations have been carried out for the case 
of universal soft terms. This is not only the most simple possibility in the framework of 
SUGRA, but is also allowed in the context of superstring models. This is e.g. the case 
of the dilaton-dominated SUSY-breaking scenario [E2f or weakly and strongly coupled 
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Figure 9: The same as in Fig. |3|b but for two values of the initial scale, Mj = 10 16 GeV 
and 10 11 GeV, and for tan /3 — 3 and 10. 



heterotic models with one Kahler modulus ||43|| . In this sense the analysis of neutralino 



nucleon cross sections of those models is included in the above analyses. This is exactly 
true for Mj = Mqut- For intermediate unification scale, note that we are assuming 
gaugino mass universality at the high energy scale, although in this scenario gauge 
couplings do not unify. This situation is in principle possible in generic supersymmetric 
models, however it is not so natural in supersymmetric models from supergravity where 
gaugino masses and gauge couplings are related through the gauge kinetic function. 
Since an explicit string construction with nonuniversal gauge couplings and gaugino 
masses will be analyzed in detail below, we have chosen to simplify the discussion here 
assuming gaugino mass universality. 

Obviously, following the discussion of Subsection 3.3, non universality of the soft 
terms in addition to intermediate scales may introduce more flexibility in the compu- 
tation. In particular, decreasing in order to obtain regions in the parameter space 
giving rise to cross sections compatible with the sensitivity of current detectors, may 
be easier. 

D-brane scenarios 

D-brane constructions are explicit scenarios where both situations mentioned above, 
non-universality and intermediate scales, may occur. The first attempts to study dark 
matter within these constructions were carried out in scenarios with the unification 
scale Mqut ~ 10 16 GeV as the initial scale E3, [TO, and dilaton-dominated SUSY- 



breaking scenarios with an intermediate scale as the initial scale [1131. However, the 
important issue of the D-brane origin of the U(l)y gauge group as a combination of 
other U(l)'s and its influence on the matter distribution in these scenarios was not 
included in the above analyses. When this is taken into account, interesting results 
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are obtained In particular, scenarios with the gauge group and particle content of 
the SUSY standard model lead naturally to intermediate values for the string scale, in 
order to reproduce the value of gauge couplings deduced from experiments. In addition, 
the soft terms turn out to be generically non universal. Due to these results, large cross 
sections in the small tan {3 regime can be obtained. 



Let us consider for example a type I string scenario where the gauge group 
£7(3) x £7(2) x £7(1), giving rise to SU(3) x SU{2) x £7(1) 3 , arises from three different 
types of D-branes, and therefore the gauge couplings are non-universal as in Fig. ^a g . 
Other examples with the standard model gauge group embedded in D-branes in a 
different way can be found in ref.fnj. Here £7(l)y is a linear combination of the 
three £7(1) gauge groups arising from £7(3), £7(2) and U(l) within the three different 
D-branes. This implies 



+ 



1 



+ 



ay(Mj) ai(Mj) a 2 (Mj) 3a 3 (Mj) 



(13) 



where a>k correspond to the gauge couplings of the U(k) branes. As shown in ref.jnj, 
ai(Mj) = 0.1 leads to the string scale Mj = 10 12 GeV. On the other hand, the extra 
£7(l)'s are anomalous and therefore the associated gauge bosons have masses of the 
order of Mj. 

The analysis of the soft terms has been done under the assumption that only the 
dilaton (S) and moduli (Tj) fields contribute to SUSY breaking and it has been found 
that these soft terms are generically non-universal. Using the standard parameteriza- 
tion ra 



F s = s/3{S + S*)m 3/2 sm9 , 
F* = V3(Ti + T*)m 3/2 cos 9 Bi 



(14) 



where i = 1,2,3 labels the three complex compact dimensions, and the angle 9 and 
the Qi with J2i |©i| 2 = 1> just parameterize the direction of the goldstino in the S, 
Ti field space, one is able to obtain the following soft terms [13]. The gaugino masses 
associated to the three gauge groups of the standard model are given by 



M 3 
M 2 



M 



Y 



y/3m 3 / 2 sm9 , 
V3m 3 /2 0i cos 9 

V3m 3/2 a Y (Mj) 



'2 3 cos# 0icos# 2 sin# 



Qti (Mj) a 2 (M 7 ) 3a 3 (Mj 



(15) 



The soft scalar masses of the three families are given by 



m l/2 



g To be precise, the running of the U{l)y gauge coupling is not exactly like in the figure since, due 
to the D-brane origin of the 17(1) gauge groups, relation ( fl3| ) must be fulfilled. See in this respect 
Fig. 2 in ref.gl. 
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where e.g. w# denotes the three family squarks ur, cr, £r. Finally the trilinear param- 
eters of the three families are 



A- 



—m 3 / 2 [(0 2 - 0i - 3 ) cos# - sin9 



2 
. 



m 3/2 [(03 — ©i — ©2) cos# — sin9 



;i7) 



Although these formulas for the soft terms imply that one has in principle five free 
parameters, m 3 / 2 , 9 and 0j with i = 1,2,3, due to relation J2i |0i| 2 — 1 om y f° ur of 
them are independent. In the analysis the parameters 9 and 0j are varied in the whole 
allowed range, < 9 < 2n, — 1 < 0j < 1. For the gravitino mass, 7723/2 < 300 GeV 
is taken. Concerning Yukawa couplings, their values are fixed imposing the correct 
fermion mass spectrum at low energies, i.e., one is assuming that Yukawa structures of 
D-brane scenarios give rise to those values. 

Fig. [XOj displays a scatter plot of function of the neutralino mass m^o for 

a scanning of the parameter space discussed above. Two different values of tan f3, 10 
and 15, are shown. LEP and Tevatron bounds on SUSY masses are included as in the 
previous sections. They forbid e.g. values of 7713/2 smaller than 170 GeV. Although 
bounds coming from CLEO b — > 37 branching ratio measurements are not included in 
the figures, one can check explicitly that their qualitative patterns are not modified. 
It is worth noticing that for tan/3 = 10 there are regions of the parameter space 
consistent with DAMA limits. In fact, one can check that tan/3 > 5 is enough to 
obtain compatibility with DAMA. Since the larger tan f3 is, the larger the cross section 
becomes, for tan /3 = 15 these regions increase. 

Let us recall that both plots in the figure are obtained taking m 3 /2 < 300 GeV, 
which corresponds to squark masses smaller than 500 GeV at low energies. Larger 
values of 7n 3 / 2 will always produce cross sections below DAMA limits. In particular, 
the right hand side and bottom of the plots will also be filled with points. Cross sections 
below projected GENIUS limits will be possible for both figures. On the other hand, 
it is worth mentioning that the isolated points in the plots with, in general, very large 
values of the cross section correspond to values of the lightest stop mass extremely 
close to the mass of the LSP, in particular (mf — mg>)/m% < 0.01. 
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Figure 10: The same as in Fig. |3|b but for the D-brane scenario with the string scale 
Mj = 10 12 GeV discussed in the text, and for tan/3 = 10 and 15. 



Finally, let us mention that scenarios where all gauge groups of the standard model 
are embedded within the same set of D-branes, and therefore with gauge coupling 
unification, are also possible. However, unlike the previous scenario, now tan f3 > 20 is 



necessary in order to obtain regions consistent with DAMA limits |TJ]]. As discussed 
above in the context of mSUGRA this is due to the different values of the a's at the 
string scale in both types of scenarios. 



4 Relic neutralino density versus cross section 

As discussed in the Introduction, current dark matter detectors are sensitive to a 
neutralino-proton cross section around 10 -6 pb. This value is obtained taking into 
account, basically, that the density of dark matter in our Galaxy, which follows from 
the observed rotation curves, is pdm ~ 0.3 GeV/cm 3 . Thus in this work we were 
mainly interested in reviewing the possibility of obtaining such large cross sections in 
the context of SUSY scenarios. In order to compute the cross section only simple field 
theory techniques are needed, no cosmological assumptions about the early Universe 
need to be used. 

On the other hand, such cosmological assumptions indeed must be taken into ac- 
count when computing the amount of relic neutralino density arising from the above 
scenarios. Generically, one obtains |J 

n 9 oh 2 ~ — - — , (is) 

where a^S n is the cross section for annihilation of a pair of neutralinos into standard 
model particles, v is the relative velocity between the two neutralinos, and < .. > 
denotes thermal averaging. The constant C involves factors of Newton's constant, the 
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temperature of the cosmic background radiation, etc. Then one may compare this 
result with dark matter observations in the Universe. Let us then discuss briefly the 
effect of relic neutralino density bounds on cross sections. 

The most robust evidence for the existence of dark matter comes from relatively 
small scales. Lower limits inferred from the flat rotation curves of spiral Galaxies || 46 



are Qdm h > 0.01 — 0.05, where h is the reduced Hubble constant. On the opposite 



side, observations at large scales, (6 — 20) h Mpc, have provided estimates J47[ of 
^DMh 2 ~ 0.1 — 0.6, but values as low as Qdm^ 2 ~ 0.02 have also been quoted Tak- 
ing up-to-date limits on h, the baryon density from nucleosynthesis and overall matter- 
balance analysis one is able to obtain a favoured range |^|, |4£| [5(J, 0.01 < VlDAih 2 < 0.3 
(at ~ 2cr CL). Note that conservative lower limits in the small and large scales are of 
the same order of magnitude. 

As is well known, for cr|o n of the order of a weak-process cross section, Q^o obtained 
from eq.(^) is within the favoured range discussed above ||. This is precisely the 
generic case when the lightest neutralino is mainly bino. Then, the neutralino-nucleus 
cross section is of the order of 1 pb, i.e. <y^o_ p « 10~ 8 pb, and therefore it is natural 
to obtain that neutralinos annihilate with very roughly the weak interaction strength. 
These cross sections were discussed in Section 2 where low and moderate values of tan j3 
within mSUGRA were considered. In fact, for these cross sections, there is always a 
set of parameters which yield 0.1 < Q^oh 2 < 0.3. This analysis, including a complete 
treatment of coannihilations was carried out in refs.[l6|, |51[| . 

On the other hand, in this review we were interested in larger neutralino-nucleon 
cross sections in order to be in the range of sensitivity of current dark matter detectors. 
It is then expected that such high neutralino-proton cross sections <Jy$- p « 10~ 6 pb, as 
those presented in Section 3, will correspond to relatively low relic neutralino densities. 
This is in general the situation. There is always a set of parameters for which the value 
of the relic density is still inside the ranges we considered above when discussing the 
observational bounds, but it is generically close to the conservative lower bound. The 
analysis of the relic neutralino density for scenarios with large tan (3 and non-universal 
scalar masses was carried out in refs.fT], ||, |, |52| . Scenarios with non- universal gaugino 



masses were studied in refs. |26|, 2%. An analysis of the relic density for focus point 



supersymmetry can be found in ref.||TT|. Discussions for scenarios with intermediate 



unification scales can be found in refs. [|P2] |-||14|p|. 

Finally, in case of preferring the stronger lower bound Q^Mh 2 > 0.1, let us mention 
the possibility that not all the dark matter in our Galaxy are neutralinos. Then Q^o < 
Qdm, and therefore Q^o < 0.1 is possible. However, due to the corresponding reduction 
in the density of neutralinos in the Galactic halo, the neutralino-proton cross section 



should be increased in order to maintain the experimental detection rates [16, 55, 56 



h Let us remark, however, that scenarios with intermediate scales might give rise to cosmological 
results different from the usual ones summarized in eq. ( |i~§| ) |53| . 
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5 Final comments and outlook 



There is overwhelming evidence that most of the matter in the Universe is dark matter. 
In the present paper we have reviewed the direct detection of supersymmetric dark 
matter in the light of recent experimental efforts. In particular, DAMA collaboration 
using a Nal detector has reported recently 0] data favouring the existence of a WIMP 
signal in their search for annual modulation. They require a large cross section of the 
order of 10~ 6 pb. We have observed that there are regions in the parameter space of 
SUGRA scenarios where such a value can be obtained, although it is fair to say 



that smaller values can also be obtained and even more easily. The latter result may 
be important since CDMS collaboration using a germanium detector has reported a 
null result for part of the region explored by DAMA. Clearly, more sensitive detectors 
producing further data are needed to solve this contradiction. Fortunately, many dark 
matter detectors are being projected. This is the case e.g. of DAMA 250 kg. and 



CDMS Soudan, but particularly interesting is the projected GENIUS |20| detector 
where values of the cross section as low as 10~ 9 pb will be accesible. 

In summary, underground physics as the one discussed here in order to detect dark 
matter is crucial. Even if neutralinos are discovered at future particle accelerators 
such as LHC, only their direct detection due to their presence in our galactic halo will 
confirm that they are the sought-after dark matter of the Universe. 

Acknowledgments 

We thank E. Gabrielli and E. Torrente-Lujan as co-authors of some of the works re- 
ported in this review. D.G. Cerdeno acknowledges the financial support of the Co- 
munidad de Madrid through a FPI grant. The work of S. Khalil was supported by 
PPARC. The work of C. Munoz was supported in part by the Ministerio de Ciencia y 
Tecnologfa, and the European Union under contract HPRN-CT-2000-00148. 

References 



[1] For a recent review, see, e.g., Y. Sofue and V. Rubin, |astro-ph /00 1 0594] , and 
references therein. 

[2] For a review, see, e.g., B. Sadoulet, Rev. Mod. Phys. 71 (1999) 197, and references 
therein. 



[3] For a review, see, e.g., M. Drees, [hep-ph / 980423 1| , and references therein. 

[4] R. Bernabei et al., Phys. Lett. B480 (2000) 23. 

[5] R. Abusaidi et al., Phys. Rev. Lett. 84 (2000) 5699. 

[6] For a review, see, e.g., G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 
267 (1996) 195, and references therein. 

[7] A. Bottino, F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D59 (1999) 095004. 



20 



E. Accomando, R. Arnowitt, B. Dutta and Y. Santoso, Nucl. Phys. B585 (2000) 
124. 



[9] M.E. Gomez and J.D. Vergados, |hep-ph/0012"020 



[10] A. Corsetti and P. Nath, |hep-ph/0003T86 



[11] J.L. Feng, K.T. Matchev and F. Wilczek, Phys. Lett. B482 (2000) 388. 

[12] E. Gabrielli, S. Khalil, C. Munoz and E. Torrente-Lujan, Phys. Rev. D63 (2001) 
025008. 

[13] D. Bailin, G.V. Kraniotis and A. Love, Phys. Lett. B491 (2000) 161. 

[14] D.G. Cerdeno, E. Gabrielli, S. Khalil, C. Munoz and E. Torrente-Lujan, Nucl. 
Phys. B603 (2001) 231. 

[15] see e.g.: A. Kudo and M. Yamaguchi, [hep-ph/0103272| , and references therein. 

[16] J. Ellis, A. Ferstl and K. Olive, Phys. Lett. B481 (2000) 304. 

[17] L.E. Ibanez and C. Lopez, Nucl. Phys. B233 (1984) 511; L.E. Ibanez, C. Lopez 
and C. Munoz, Nucl. Phys. B256 (1985) 218. 



[18] For recent reviews, see, e.g., A. Czarnecki and W.J. Marciano, [hep-ph /0 102 122| ; 
S. Narison, |hep-ph/0103T99| . 



[19] DELPHI collaboration, [hep-ex/0104020| , 0104021. 



[20] H.V. Klapdor-Kleingrothaus, L. Baudis, G. Heusser, B. Majorovits and H. Paes, 
[hep-ph/9910205 . 



[21] J.L. Feng, K.T. Matchev and T. Moroi, Phys. Rev. Lett. 84 (2000) 2322, Phys. 
Rev. D61 (2000) 75005. 



[22] For a review, see, e.g., A. Brignole, L.E. Ibanez and C. Munoz, |hep-ph / 9707209| , 
and references therein. 

[23] M. Olechowski, S. Pokorski, Phys. Lett. B344 (1995) 201. 

[24] V. Berezinsky, A. Bottino, J. Ellis, N.Fornengo, G. Mignola and S. Scopel, As- 
tropart. Phys. 5 (1996) 1, ibid, 5 (1996) 333. 

[25] R. Arnowitt and P. Nath, Phys. Rev. D56 (1997) 2820. 

[26] S. Mizuta, D. Ng and M. Yamaguchi, Phys. Lett. B300 (1993) 96. 

[27] A. Birkedal-Hansen and B. Nelson, Phys. Rev. D64 (2001) 015008. 

[28] T. Moroi and L. Randall, Nucl. Phys. B570 (2000) 455. 



21 



[29 
[30 
[31 

[32 

[33 
[34 
[35 
[36 
[37 
[38 

[39 

[40 
[41 
[42 



[43 
[44; 
[45 
[46 
[47 
[48 
[49 



B. Murakami and J.D. Wells, |hep-pri/00ll08^ 
J. Lykken, Phys. Rev. D54 (1996) 3693. 

N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. Lett. B249 (1998) 262; 
I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 
(1998) 263; I. Antoniadis and C. Bachas, Phys. Lett. B450 (1999) 83. 

G. Shiu and S.-H.H. Tye, Phys. Rev. D58 (1998) 106007; Z. Kakushadze and 
S.-H.H. Tye, Phys. Rev. D58 (1998) 126001. 

K. Benakli, Phys. Rev. D60 (1999) 104002. 

C. Burgess, L.E. Ibanez and F. Quevedo, Phys. Lett. B447 (1999) 257. 
I. Antoniadis and B. Pioline, Nucl. Phys. B550 (1999) 41. 

K. Benakli and Y. Oz, Phys. Lett. B472 (2000) 83; A. Gregori, |hep-th/0005T9l . 
For a review, see, e.g., C. Munoz, [hep-ph /9709329 , and references therein. 



S.A. Abel, B.C. Allanach, F. Quevedo, L.E. Ibanez and M. Klein, JHEP 0012 
(2000) 026. 

V.A. Kuzmin and V.A. Rubakov, Phys. Atom. Nucl. 61 (1998) 1028; M. Birkel 
and S. Sarkar, Astropart. Phys. 9 (1998) 297. 

N. Kaloper and A. Linde, Phys. Rev. D59 (1999) 101303. 

L.E. Ibanez, C. Munoz and S. Rigolin, Nucl. Phys. B553 (1999) 43. 



L.E. Ibanez, |hep-ph / 9905349| ; I. Antoniadis, C. Bachas and E. Dudas, Nucl. Phys. 
B560 (1999) 93; N. Arkani-Hamed, S. Dimopoulos and J. March-Russell, [hep-| 
th/9908l46| . 



For a review, see, e.g., C. Munoz, [hep-th/9906152 , and references therein. 
S. Khalil, Phys. Lett. B484 (2000) 98. 



R. Arnowitt, B. Dutta and Y. Santoso, [hep-ph/0005154 
P. Salucci and M. Persic, |astro-ph/ 9703021 , 9903432. 



see e.g.: W.L. Freedman, |astro-ph/ 9905222] , and references therein. 



N. Kaiser, |astro-ph/ 9809341 



J.R. Primack, |astro-ph / 00071 S7\ J.R. Primack and M.A.K. Gross, |astro 



ph/0007165 



[50] K. Griest and M. Kamionkowski, Phys. Rep. D64 333-334 (2000) 167. 



22 



[51] J. Ellis, T. Falk, G. Ganis, K.A. Olive and M. Srednicki, |hep-ph/0l02Qgg 



[52] R. Arnowitt, B. Dutta and Y. Santoso, |hep-ph/0102T81 
[53] S. Khalil et al., in preparation. 

[54] G.F. Giudice, E.W. Kolb and A. Riotto, Phys. Rev. D64 (2001) 023508. 
[55] J. Ellis, A. Ferstl and K. Olive, |hep-ph/0007Tl3| . 
[56] Y.G. Kim and M. Nojiri, |hep-ph/ 01 04258 . 



23 



